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1.  INTRODUCTION 


Over  the  past  several  years  our  research  efforts  under  this  contract  have 
been  focused  upon  developing  techniques  to  utilize  satellite  data  assimilation 
in  order  to  improve  the  forecasts  made  by  global  numerical  weather  prediction 
models.  The  accuracy  of  forecasts  made  by  forecast  models  is  dependent 
upon  how  realistically  the  model  represents  the  actual  processes  of  the 
atmosphere  and  by  the  ability  to  provide  the  model  with  initial  conditions 
which  reflect  the  true  state  of  the  atmosphere.  Data  assimilation  research 
has  concentrated  on  the  second  of  these  two  problems.  Before  the  develop¬ 
ment  of  satellite-borne  measurement  systems  in  the  1960's,  the  volume  of 
asynoptic  data  available  was  limited  and  data  assimilation  schemes  were 
designed  primarily  to  handle  the  conventional  surface  and  upper-air  obser¬ 
vations.  Unlike  the  conventional  data,  the  vast  amount  of  satellite  observa¬ 
tions  consists  of  measurements  of  meteorological  variables  (primarily  temper¬ 
ature)  which  are  widely  distributed  in  both  space  and  time.  Beginning  with 
a  paper  by  Charney  et  al_^  (1969),  a  great  deal  of  data  assimilation  research 
has  been  conducted  to  determine  ways  to  use  this  new  data  source.  In  this 
report  we  will  use  the  terminology  conventions  introduced  in  the  review 
papers  of  McPherson  (1975)  and  Bengtsson  (1975). 

Dynamic  assimilation  describes  the  useof  a  numerical  prediction  model 
as  an  integrator  of  observations  distributed  in  space  and  time.  We  refer 
to  the  process  of  interpolating  an  observation  to  a  number  of  surrounding 
model  grid  points  using  an  objective  analysis  scheme  as  indirect  insertion. 
Insertion  of  data  into  a  numerical  model  which  disagrees  with  the  current 
model  solution  will  result  in  the  generation  of  gravity  wave  noise  within  the 
model.  On  the  other  hand  insertion  of  data  which  is  in  perfect  agreement  with 


the  model  solution  adds,  no  information.  There  are  two  ways  to  attack  the 
data  assimilation  problem;  use  an  insertion  method  which  maximizes  infor¬ 
mation  and  minimizes  resulting  model  noise  or  control  the  model  noise  after 
insertion.  Ghi  1  et  al^  (.1979.)  found  in  a  global  assimilation  experiment 
using  real  data  that  forecast  improvement  was  directly  related  to  the 
sophistication  of  the  continuous  indirect  insertion  scheme  used.  On  the 
other  hand,  in  our  research  we  have  utilized  techniques  (Noise  Freezing 
Methods  I  and  II)  to  control  the  noise  caused  by  the  insertion.  In  the 
research  described  in  Sasaki  and  Goerss  (.1982)  we  used  a  dry  global  primi¬ 
tive  equations  model  with  3  vertical  levels  and  7.5°  resolution  and  a  rela¬ 
tively  simple  indirect  insertion  scheme.  Continuous  assimilation  of  sate¬ 
llite  temperatures  was  performed  using  Noise  Freezing  Methods  I  and  II  and 
the  results  were  compared  with  a  scheme  designed  to  simulate  most  operational 
forecast  centers.  Significant  forecast  improvement  was  found  when  the  noise 
control  techniques  were  utilized. 

>In  this  report  we  describe  the  culmination  of  our  global  satellite  data 
assimilation  research  efforts.  A  continuous  temperature  data  assimilation 
scheme  was  designed  to  be  implemented  using  the  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS).  In  the  next  section  the  details 
of  this  scheme  are  described.  Two  experiments  were  run  using  the  operational 
satellite  temperature  data  available  to  Fleet  Numerical  Oceanographic  Center 
(FNOC).  Each  experiment  consisted  of  an  assimilation  and  a  control  run.  The 
results  of  these  experiments  are  detailed  in  Section  3.  We  conclude  our 
report  by  summarizing  our  research  and  discussing  its  future  applications 
to  the  improvement  of  NOGAPS  forecasts. 
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2.  CONTINUOUS  DATA  ASSIMILATION  TECHNIQUES 


The  key  to  successful  data  assimilation  is  to  introduce  information 
to  the  forecast  model  without  "shocking"  the  model  so  much  that  the  benefits 
of  the  added  information  are  overcome  by  the  noise  induced  by  the  insertion. 
As  described  by  Miyakoda  et  aK.  (1978),  when  asynoptic  data  are  inserted 
into  an  atmospheric  model,  they  create  discontinuity  with  the  current  model 
solution.  The  dynamical  characteristics  of  the  fluid  enable  it  to  remove 
this  discontinuity  in  such  a  way  as  to  maintain  a  state  of  approximate  geo- 
strophic  balance.  This  geostrophic  adjustment  process  results  in  the  gener¬ 
ation  and  dispersal  of  inertia  gravity  waves.  Talagrand  (1972)  points  out 
that  the  effectiveness  of  an  assimilation  scheme  is  dependent  upon  not  only 
the  amount  of  information  introduced  into  the  model  but  also  upon  how  well 
the  high  frequency  noise  caused  by  this  introduction  can  be  damped. 

The  conclusion  of  two  studies  using  simulated  asynoptic  data  (Bengtsson 
and  Gustavsson,  1971;  Miyakoda  and  Talagrand,  1971)  is  that  the  spreading  of 
information  contained  in  an  observation  to  several  surrounding  grid  points 
using  some  sort  of  objective  analysis  scheme  increases  the  effectiveness 
of  assimilation.  In  an  extensive  experiment  using  real  data  and  a  number 
of  indirect  insertion  schemes,  Ghi 1  et  £1^  (1979)  found  that  the  success  of 
assimilation  (based  on  forecast  improvement)  was  directly  related  to  the 
level  of  sophistication  of  the  insertion  scheme  used.  The  reason  for  this 
is  that  the  asynoptic  data  is  put  into  a  form  by  the  objective  analysis 
scheme  which  makes  it  more  "palatable"  to  the  forecast  model.  The  effec¬ 
tiveness  of  this  transformation  in  terms  of  forecast  improvement  is  depen¬ 
dent  upon  how  much  of  the  information  contained  in  the  asynoptic  data  is 
retained.  In  this  study  we  have  devised  a  two-step  objective  analysis 
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procedure  to  attack  this  problem.  First,  we  utilize  a  variational  adjust¬ 
ment  scheme  which  we  call  vertical  shape  matching  to  modify  each  satellite- 
derived  temperature  sounding  so  that  its  information  content  is  disturbed 
as  little  as  possible  but  its  form  is  changed  so  that  it  can  be  more  readily 
accepted  into  the  NOGAPS  model.  Then,  the  soundinqs  are  assimilated  each 
time-step  into  the  NOGAPS  model  using  an  objective  analysis  scheme  patterned 
after  that  described  by  Barker  (1982).  The  details  of  these  procedures  will 
be  given  later  in  this  section. 

Another  aspect  of  data  assimilation  has  been  explored  by  Sasaki  and 
Goerss  (1982).  In  their  study  of  Noise  Freezing  Methods  they  found  that 
improved  forecasts  could  be  obtained  utilizing  continuous  data  assimilation 
and  by  modifying  the  numerical  forecast  model  so  that  the  noise  induced  by 
data  insertion  is  controlled.  Their  work  pointed  out  the  need  to  examine 
closely  the  effects  of  the  data  insertion  upon  the  model  used  in  the  assimi¬ 
lation  process.  In  this  section  we  will  also  discuss  this  aspect  of  the 
problem  with  respect  to  the  NOGAPS  model. 

Vertical  Shape  Matching 

The  vertical  resolution  of  satellite  soundings  is  such  that  smoothing 
results  for  features  such  as  frontal  inversions  and  the  tropopause  which 
are  better  defined  by  forecast  models.  An  example  of  this  problem 
is  illustrated  in  Fig.  1  where  we  have  pictured  a  model  sounding  with  two 
inversions  and  a  satellite  sounding  in  which  these  features  are  typically 
smoothed.  Before  the  observed  satellite-derived  temperatures  are  inserted 
into  the  forecast  model,  a  variational  adjustment  is  made  to  each  satellite 
sounding.  This  adjustment  is  accomplished  by  minimizing  the  functional: 
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where  the  vertical  coordinate  is  n  =  In  (1000/p).  The  mean  layer  temperatures 
implied  by  the  satellite  sounding  are  given  by  T  while  the  model -derived 
temperatures  at  the  location  of  the  satellite  sounding  are  given  by  T.  The 
resulting  adjusted  satellite  sounding,  given  by  T,  is  such  that  the  vertical 
structure  of  the  model  is  preserved  along  with  the  mean  layer  temperatures 
implied  by  the  satellite  sounding.  The  results  of  applying  this  technique 
are  illustrated  in  Fig.  1.  Thus,  by  applying  this  scheme  we  produce  sate¬ 
llite  soundings  which  have  been  modified  in  a  way  so  that  the  information 
content  that  we  can  expect  from  such  a  sounding  is  preserved  and  yet  at  the 
same  time  the  sounding  has  been  made  more  consistent  with  the  model  into 
which  it  will  be  inserted. 


Indirect  Insertion  Scheme 

The  indirect  insertion  scheme  used  in  this  study  is  an  objective  analysis 
procedure  patterned  after  that  described  by  Barker  (1982).  For  a  particular 
time-step  the  model  temperature  fields  are  modified  by  performing  a  two-pass 
Barnes  analysis  using  all  satellite  soundings  to  be  applied  at  that  time. 

Let  K  represent  the  number  of  asynoptic  observations  taken  during  a  block  of 
time  within  radius  R  of  a  model  grid  point.  For  this  study  R  was  chosen  to 
be  800  km.  For  each  assimilation  time-step,  a  two-pass  Barnes  analysis  is 
performed.  On  the  first  pass  the  forecast  temperature  T^,  at  each  grid  point, 
is  replaced  by 

K  K 

T  =  Tf  +  [  Z  W.(T.  -  T-. )/(W  +  l  Wj], 
r  i=i  1  1  ri  o  1 
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where  T..  is  the  observed  temperature  at  the  pressure  level  of  the  grid 
point  and  T^.  is  the  model  forecast  temperature  interpolated  to  the  loca¬ 
tion  and  pressure  level  of  the  observation.  The  weight  associated  with  the 
observation  is 

Wi  =  exp  [-0.2/(Y82)], 

where  is  the  distance  from  the  grid  point  to  the  location  of  the  ith 

observation  and  8  is  chosen  to  be  450  km.  The  value  W  is  a  "background" 

o 

weight  which  results  in  the  value  of  the  forecast  field  being  given  weight 
as  if  it  were  an  observation  450  km  from  the  grid  point.  This  prevents  a 
"noisy"  analysis  at  grid  points  on  the  fringes  of  areas  where  satellite 
observations  are  taken.  The  process  is  repeated  on  the  second  pass  with 
the  value  of  y  reduced  from  1.0  to  0.3  and  with  the  forecast  values  being 
replaced  by  the  results  of  the  first  pass. 

Continuous  Data  Assimilation  into  NOGAPS 


The  main  forecast  tool  of  Fleet  Numerical  Oceanographic  Center  is  the 
Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS).  NOGAPS 
consists  of  objective  analysis,  initialization,  and  forecast  model.  The 
NOGAPS  model  is  with  minor  modifications  a  version  of  the  UCLA  general  cir¬ 
culation  model  described  by  Arakawa  and  Lamb  (1977).  Specific  details  about 
the  NOGAPS  model  are  given  by  Rosmond  (1981).  Details  on  the  objective 
analysis  method  and  initialization  procedure  are  given  by  Barker  (1981). 

In  this  study  we  utilized  a  version  of  the  NOGAPS  model  for  which  the 
time-step  was  240  seconds  and  one  "cycle"  of  integration  consisted  of  a 
Matsuno  time-step  and  four  leap-frog  time-steps.  The  last  of  the  leap-frog 
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time-steps  included  most  of  the  physical  parameter izations  of  diabatic 
processes.  The  Matsuno  time-step  has  a  damping  effect  upon  the  "noise" 
induced  within  the  forecast  model  by  the  full  physics  computations  performed 
every  five  time-steps.  In  Fig.  2  we  see  the  sawtooth  pattern  this  cycle 
produces  in  the  global  root  mean  square  pressure  tendency  for  a  six-hour 
forecast  beginning  at  06Z,  27  March  1983.  Over  the  six-hour  period  the 
level  of  this  quantity  is  virtually  constant  at  roughly  0.0006  mb/sec. 

The  satellite  data  used  in  this  study  were  those  routinely  available 
on  an  operational  basis  at  FNOC.  The  data  we  used  consisted  of  observations 
of  geopotential  height  at  the  ten  standard  pressure  levels  between  1000  mb 
and  100  mb,  inclusive.  Using  the  hypsometric  equation  these  observations 
were  converted  to  a  mean  layer  temperature  which  was  then  assumed  to  repre¬ 
sent  the  atmospheric  temperature  at  the  pressure  level  at  the  middle  of  the 
layer.  Thus,  for  each  satellite  sounding  we  obtained  nine  temperatures  at 
pressure  levels  ranging  from  922  mb  to  122  mb,  inclusive.  We  then  blocked 
the  satellite  data  into  one-hour  blocks  centered  upon  the  hour  of  observation. 
Typically  each  block  consisted  of  around  100  observations  and  global  cov¬ 
erage  was  effected  by  about  twelve  hours  worth  of  data.  The  vertical  shape 
matching  discussed  previously  in  this  section  was  performed  for  each  block 
of  satellite  data  immediately  before  that  block  was  introduced  into  the 
model.  Thus,  for  a  block  of  satellite  data  centered  on  09Z  the  model  fields 
valid  at  0830Z  were  used  along  with  the  satellite  observations  themselves 
in  order  to  modify  the  soundings  which  would  be  inserted  into  the  model 
between  0830Z  and  0930Z.  The  data  insertion  was  performed  during  each  of 
the  three  interior  leap-frog  time-steps  of  each  N0GAPS  cycle. 

As  mentioned  previously  Sasaki  and  Goerss  (1982)  found  that  forecast 


improvement  from  data  ass/imlation  was  dependent  upon  how  well  the  numer¬ 
ical  forecast  model  controlled  the  noise  induced  by  data  insertion.  For  the 
forecast  model  they  used  in  that  study,  they  found  that  data  assimilation 
without  use  of  one  of  the  Noise  Freezing  Methods  resulted  in  a  gradual 
increase  in  the  global  root  mean  square  pressure  tendency  over  the  forecast 
period.  The  use  of  Noise  Freezing  Method  I,  in  which  a  constant  was  sub¬ 
tracted  from  all  values  of  surface  pressure  carried  within  the  model,  resulted 
in  the  same  gradual  increase  in  this  quantity  but  at  a  substantially  reduced 
level.  The  results  of  applying  Noise  Freezing  Method  II,  in  which  periodically 
during  the  assimilation  process  the  mass  divergence  fields  were  smoothed  and 
the  model  winds  adjusted  to  agree  with  the  smoothed  mass  divergence,  showed 
a  sawtooth  pattern  in  which  the  level  of  the  global  rms  pressure  tendency 
is  constant  but  higher  than  the  level  of  that  quantity  for  a  non-assimilation 
forecast.  In  Fig.  2,  we  can  see  that  the  results  obtained  for  a  six-hour 
assimilation  forecast  beginning  at  06Z,  27  March  1983  are  similar  to  those 
obtained  by  Sasaki  and  Goerss  (.1982)  using  Noise  Freezing  Method  II.  The 
Matsuno  time-step  every  cycle  provides  the  noise  control  necessary  for  the 
data  assimilation  forecast. 

One  final  aspect  of  data  assimilation  will  now  be  discussed.  There  are 
two  extremes  which  can  be  considered  when  performing  data  assimilation.  If 
the  information  to  be  added  is  directly  inserted  into  the  forecast  model 
with  no  modification  whatsoever,  the  resulting  noise  induced  within  the 
model  will  be  maximal.  On  the  other  hand,  if  the  information  is  modified 
so  much  that  it  agrees  perfectly  with  the  model  solution  at  time  of  insertion, 
then  no  noise  will  be  induced  but  no  information  will  be  added  to  the  model. 
Recall  that  our  assimilation  procedure  calls  for  the  satellite  data  to  be 
broken  up  into  one-hour  blocks  centered  on  the  hour  of  observation.  During 


each  of  these  hlocks,  the  data  Is  inserted  into  the  model  on  each  of  the 
three  Interior  leap-frog  ti.me-s.teps  of  each,  model  cycle..  In  Fig.  2  the 
astericks  represent  the  beginning  of  a  new  block  of  satellite  information. 

At  the  top  of  the  figure  Is  plotted  the  root  mean  square  of  the  temperature 
differences  between  the  model  solution  for  that  time-step  at  the  location 
of  the  observation  and  the  value  of  the  observation  itself.  This  is  for 
the  satellite  observations  for  temperature  at  about  592  mb  and  these  results 
are  typical  of  those  found  at  all  levels  except  for  the  uppermost  level 
(about  122  mb)  where  the  rms  was  typically  three  to  four  times  greater.  The 
important  things  to  note  from  this  plot  are  that  the  satellite  information 
is  different  from  the  model  solution  (as  indicated  by  the  high  rms  values 
for  the  first  two  time-steps)  and  that  the  information  is  being  successfully 
introduced  into  the  model  (as  indicated  by  the  low  rms  values  for  the  suc¬ 
ceeding  time-steps).  Thus,  we  conclude  that  our  data  assimilation  proce¬ 
dures  are  sound  and  proceed  with  the  experiment  itself. 


3.  EXPERIMENTAL  RESULTS 


Two  continuous  data  assimilation  experiments  were  conducted,  the  first 
between  May  12-18,  1983  and  the  second  between  August  22-27,  1983.  Each 
experiment  consisted  of  an  assimilation  run  and  a  control  run  using  NOGAPS 
and  the  satellite  data  operationally  available  to  FNOC.  A  six-hour  update 
cycle  was  utilized  in  these  experiments.  Each  experiment  consisted  of  a 
start-up  period  climaxed  by  72-hour  forecasts  made  from  the  global  fields 
produced  by  \.he  assimilation  and  control  runs.  In  this  section  we  will 
describe  the  differences  between  an  assimilation  and  a  control  run,  outline 
how  each  experiment  was  conducted,  and  evaluate  the  results  of  the  forecasts 
made  at  the  end  of  each  experiment. 

Both  assimilation  and  control  runs  used  NOGAPS  as  described  by  Rosmond 
(1981)  and  Barker  (1981).  The  version  of  the  forecast  model  used  here  pos¬ 
sessed  six  vertical  layers,  2.4°  x  3.0°  horizontal  resolution,  and  a  4  minute 
time  step.  At  six-hour  intervals  objective  analysis  and  initialization  were 
performed  using  all  data  operationally  available  at  that  time  including 
satellite  observations.  The  only  difference  between  the  assimilation  run 
and  the  control  run  was  the  first-guess  fields  used  by  the  objective  analysis 
procedure  each  six  hours.  For  the  control  run  these  fields  were  simply  the 
six-hour  model  forecast  valid  at  that  time.  For  the  assimilation  run  these 
fields  were  the  result  of  re-running  that  six-hour  model  forecast  while  the 
satellite  observations  collected  over  that  time  period  were  continuously 
assimilated  into  the  model.  Thus,  the  additional  expense  of  utilizing  a 
continuous  data  assimilation  scheme  in  an  operational  setting  would  be  that 
of  an  extra  six-hour  forecast  each  update  cycle. 

Next  we  will  describe  how  each  experiment  was  conducted.  Beginning  with 
the  same  set  of  initial  fields  from  18Z,  12  May  1983  both  assimilation  and 


conventional  forecasts,  were  made  producing  fields  valid  at  00Z,  13  May  1983. 
The  satellite  observations  taken  between  18Z  and  00Z  were  continuously  assi¬ 
milated  into  the  forecast  model  producing  the  assimilation  forecast.  Using 
the  00Z  data  and  the  respective  first-guess  fields,  objective  analysis  and 
initialization  were  performed  for  both  the  assimilation  and  control  runs. 

This  process  was  repeated  each  six  hours  until  objectively  analyzed  and  ini¬ 
tialized  fields  were  produced  for  00Z,  15  May  1983.  Over  the  course  of  this 
two-day  period  the  fields  (especially  temperature)  for  the  assimilation  run 
evolved  quite  differently  than  those  for  the  control  run.  Finally,  72- 
hour  forecasts  were  made  using  the  conventional  forecast  model  using  the  00Z, 
15  May  1983  fields  from  the  assimilation  and  control  runs.  The  second  experi¬ 
ment  was  conducted  in  an  identical  fashion  beginning  with  the  same  set  of 
initial  fields  from  18Z,  21  August  1983  and  ending  with  72-hour  forecasts 
made  using  the  00Z,  24  August  1983  fields  from  the  assimilation  and  control 
runs. 

One  of  the  most  striking  results  of  these  experiments  was  the  marked 
changes  effected  in  the  temperatures  in  the  upper  levels  of  the  NOGAPS  fore¬ 
cast  model  by  continuous  temperature  data  assimilation.  During  the  course 
of  the  two-day  start-up  period  in  each  experiment  the  global  mean  temperatures 
for  the  four  lowest  model  layers  were  virtually  unaffected  by  data  assimila¬ 
tion  while  those  for  the  top  two  layers  were  significantly  altered.  This  is 
illustrated  in  Table  1  where  the  results  from  the  second  experiment  are 
tabulated.  Similar  results  were  found  from  the  first  experiment. 


TABLE  1.  Global  Mean  Temperatures  for  NOGAPS  Model  Layers 


I 


Model  Layer 


Time 

Run 

1 

2 

3 

4 

5 

6 

83082200 

A 

C 

-50.3 

-49.5 

-65.0 

-65.1 

-39.9 

-40.0 

-14.0 

-14.1 

0.6 

0.5 

12.0 

11.9 

83082212 

A 

C 

-53.9 

-49.4 

-64.0 

-64.8 

-40.1 

-39.9 

-14.0 

-14.1 

0.7 

0.4 

12.3 

11.9 

83082300 

A 

C 

-56.3 

-49.6 

-63.4 

-64.9 

-40.4 

-39.9 

-14.0 

-14.1 

0.7 

0.3 

12.3 

11.8 

83082312 

A 

C 

-59.0 

-49.5 

-62.7 

-64.8 

-40.5 

-39.8 

-13.9 

-14.1 

0.8 

0.3 

12.6 

11.9 

83082400 

A 

C 

-60,5 

-49.5 

-62.3 

-64.8 

-40.6 

-39.9 

-13.9 

-14.2 

0.8 

0.2 

12.4 

11.8 

As  can  be  seen  in  the  table  the  mean  temperatures  for  the  top  two  layers  stay 
relatively  constant  throughout  the  period  for  the  control  run  while  the  top 
layer  (^75  mb)  cools  over  10°C  and  the  second  layer  (^150  mb)  warms  almost 
3°C  for  the  assimilation  run.  More  modest  temperature  changes  are  recorded 
for  the  lower  four  layers.  The  result  is  a  more  realistic  representation  of 
the  atmospheric  temperature  profile  in  the  upper  layers  of  the  model. 

In  order  to  evaluate  the  impact  of  data  assimilation  upon  the  NOGAPS 
forecasts,  root  mean  square  errors  were  computed  for  the  differences  between 
actual  observations  and  the  value  of  the  initial  or  forecast  field  interpolated 
to  the  location  of  the  observation.  Tables  2  and  3  illustrate  the  results  of 
these  computations  for  the  first  experiment.  The  initial  height  fields  at 
various  levels  at  00Z,  15  May  1983  for  the  control  and  assimilation  runs, 
respectively,  are  featured  in  these  tables.  Globally,  we  can  see  that  the 
fields  from  the  assimilation  run  display  rms  improvement  which  increases  with 


height  and  is  most  pronounced  at  10Q  mh.  This  pattern  holds  true  in  both 
hemispheres  as  well  with  the  most  marked  improvement  displayed  at  100  mb  for 
the  Southern  Hemisphere  (over  20%  reduction  in  the  rmse).  At  850  mb  the 
assimilation  fields  globally  and  in  each  hemisphere  show  a  slight  rms  de¬ 
gradation.  The  improvements  are  also  greater  in  the  Southern  Hemisphere 
than  they  are  in  the  Northern  Hemisphere  for  the  levels  above  850. 

The  results  of  rmse  computations  for  initial  temperature  fields  are 
shown  in  Tables  4  and  5.  At  both  500  mb  and  100  mb  the  assimilation  fields 
show  large  rms  reduction  globally  and  in  each  hemisphere.  The  slight  global 
improvement  shown  at  850  mb  is  entirely  due  to  the  improvement  in  the 
Northern  Hemisphere.  Slight  degradation  is  seen  uniformly  at  300  mb.  A 
large  part  of  the  reduction  at  100  mb  is  due  to  the  reduction  in  bias  ef¬ 
fected  by  satellite  data  assimilation.  Recall  how  the  NOGAPS  model  tends 
to  be  too  warm  at  its  upper  level  without  data  assimilation.  It  should  also 
be  pointed  out  at  this  time  that  a  gross  error  check  is  performed  before  the 
rms  calculations  are  made.  At  all  levels  but  100  mb  the  number  of  rejected 
observations  is  negligible.  We  can  seen  in  Tables  4  and  5  that  nearly  150 
100  mb  observations  were  rejected  for  the  control  run  fields.  Thus,  the 
actual  rms  improvement  at  100  mb  is  even  greater  than  that  shown  in  these 
tables.  These  results  are  consistently  observed  throughout  both  experiments. 

Next  we  will  examine  the  results  for  the  24-hour  forecasts  valid  at 
00Z,  16  May  1983  displayed  in  Tables  6-9.  With  the  exception  of  the  100  mb 
level,  slight  degradation  is  shown  at  each  level  for  assimilation  height 
forecasts,  globally  and  in  each  hemisphere.  At  100  mb  significant  improve¬ 
ment  is  still  seen  (almost  20%  in  the  Southern  Hemisphere).  On  the  other 
hand  the  assimilation  temperature  forecasts  (with  the  exception  of  850  mb 
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Table  6.  Root  mean  square  error  (m)  for  24-hour  forecast  height  fields  valid 
at  00Z,  16  May  1983  for  the  control  run. 


Level  (mb] 


Global 

Obs.  RMSE 


Table  7.  Root  mean  square  error  (m)  for  24-hour  forecast  height  fields  valid 
at  00Z,  16  May  1983  for  the  assimilation  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NH 

Obs.  RMSE 


SH 

Obs.  RMSE 


Table  8.  Root  mean  square  error  (°C)  for  24-hour  forecast  temperature  fields 
valid  at  00Z,  16  May  1983  for  the  control  run. 


Level  (mb) 


Global 

Obs .  RMSE 


NH 

Obs .  RMSE 


SH 

Obs.  RMSE 


Table  9.  Root  mean  square  error  (°C)  for  24-hour  forecast  temperature  fields 
valid  at  QOZ,  16  May  1983  for  the  assimilation  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NH 

Obs.  RMSE 


SH 

Obs .  RMSE 
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and  300  mb  in  the  Southern  Hemisphere)  all  show  improvement  with  the  most 
marked  improvement  at  100  mb. 

The  results  for  the  48-hour  forecasts  valid  at  00Z,  17  May  1983  are 
shown  in  Tables  10-13.  As  for  the  24-hour  forecasts,  only  at  100  mb  is 
significant  improvement  seen  in  assimilation  height  forecasts  (roughly  10% 
in  each  hemisphere).  Only  slight  improvements  or  degradations  are  observed 
at  other  levels.  For  the  temperature  forecasts,  except  at  300  mb  where  slight 
degradations  are  seen,  improvements  are  observed  at  every  level.  Again  the 
greatest  improvement  results  at  100  mb  where  it  is  approximately  30%  glo¬ 
bally  and  even  larger  in  the  Southern  Hemisphere. 

Finally,  the  72-hour  forecast  results  are  tabulated  in  Tables  14-17. 

For  the  height  fields,  only  slight  improvements  or  degradation  are  seen  for 
the  lower  three  levels  while  significant  forecast  improvement  is  realized 
for  the  assimilation  run  at  100  mb  (about  10%  in  the  Southern  Hemisphere). 

A  similar  pattern  holds  true  for  temperature  also.  Large  forecast  improve¬ 
ment  is  observed  at  100  mb  for  the  assimilation  run  while  at  the  lower  levels 
no  significant  improvement  or  degradation  is  noted. 

In  summary,  for  the  first  forecast  experiment  we  have  seen  that,  in 
terms  of  root  mean  square  error,  the  most  marked  improvements  in  forecasts 
due  to  continuous  temperature  data  assimilation  are  seen  at  100  mb  for  both 
the  height  and  temperature  fields.  Although  the  initial  height  fields  at  the 
lower  levels  for  the  assimilation  run  showed  smaller  rmse's  than  those  for 
the  control  run,  no  significant  forecast  improvement  was  seen  at  those  levels. 
On  the  other  hand,  the  forecast  temperature  fields  at  500  mb  consistently 
showed  improvement  due  to  assimilation  while  those  at  850  mb  did  so  less 
consistently.  Little  improvement  or  degradation  was  observed  for  temperature 
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Table  12. 


Root  mean  square  error  (°C)  for  48-hour  forecast  temperature  fields 
valid  at  00Z,  17  May  1983  for  the  control  run. 


Level  (mb) 


Global 

Obs .  RMSE 


NH 

Obs.  RMSE 


SH 

Obs.  RMSE 


forecasts  at  300.  mb. 

Next  we'll  examine  the  results  for  the  second  experiment.  The  rmse's 
for  the  initial  height  and  temperature  fields  valid  at  00Z,  24  August  1983 
are  shown  in  Tables  18-21.  The  rms  improvement  for  the  heights  due  to 
assimilation  is  most  pronounced  at  500  mb  and  100  mb  with  more  modest  im¬ 
provements  at  850  mb  and  300  mb.  This  pattern  also  holds  true  for  the 
temperature  fields  with  the  largest  improvements  noted  at  100  mb.  Again, 
we'll  point  out  that  the  improvements  shown  in  these  tables  at  100  mb  are 
understated  due  to  the  number  of  rejected  observations  for  the  control  runs. 

In  Tables  22-25  are  displayed  the  results  for  the  24-hour  forecasts 
valid  at  00Z,  25  August  1983.  Small  forecast  improvements  and  degradations 
are  observed  for  the  height  fields  below  100  mb.  At  100  mb  the  assimilation 
forecast  heights  show  improvement  over  those  for  the  control  run  but  not  to 
the  degree  that  was  seen  in  the  first  experiment.  Also,  the  improvement 
seen  in  the  Northern  Hemisphere  is  larger  than  that  in  the  Southern  Hemi¬ 
sphere.  The  assimilation  temperature  forecasts  show  consistent  improvement 
with  the  most  marked  improvement  at  100  mb  for  the  Northern  Hemisphere 
(over  40%) . 

The  results  for  the  48-hour  forecasts  valid  at  00Z,  26  August  1983  are 
shown  in  Tables  26-29.  For  the  most  part  the  assimilation  height  fields  show 
slight  improvement  over  those  from  the  control  run  but  none  of  the  fields 
shows  significant  improvement.  More  consistent  improvement  is  displayed 
for  the  assimilation  forecast  temperature  fields.  Again  the  most  notable 
improvement  is  seen  at  100  mb  for  the  Northern  Hemisphere  (over  40%). 

Finally,  the  rmse  computations  for  the  72-hour  forecasts  valid  at  00Z, 

27  August  1983  are  shown  in  Tables  30-33.  Modest  improvement  for  the  height 


Table  18.  Root  mean  square  error  (m)  for  the  initial  height  fields  at  00Z, 
24  August  1983  for  the  control  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NH 

Obs .  RMSE 


SH 

Obs.  RMSE 


Table  19.  Root  mean  square  error  (m)  for  the  initial  height  fields  at  00Z, 
24  August  1983  for  the  assimilation  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NK 

Obs .  RMSE 


SH 

Obs .  RMSE 


Table  20.  Root  mean  square  error  (°C)  for  the  intial  temperature  fields  at 
00Z,  24  August  1983  for  the  control  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NH 

Obs .  RMSE 


SH 

Obs .  RMSE 


Table  21.  Root  mean  square  error  (°C)  for  the  initial  temperature  fields  at 
00Z,  24  August  1983  for  the  assimilation  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NH 

Obs .  RMSE 


SH 

Obs.  RMSE 


21.5 

28.3 

3 

40.1 

3 

65.3 

3 

recast  height 

f  i  el 

on  run. 

!0 

24 

13 

54 

8 

63 

3 

82 

10 

26 

13 

51 

8 

61 

3.32 

2.29 

3 

2.14 

3 

657 

2.96 

577 

2.84 

80 

3.74 

'  ^ -  *_• 

377 

2.13 

991 

2.04 

386 

2.34 

■ 

341 

2.24 

956 

2.04 

385 

2.68 

252 

3.67 

880 

3.40 

372 

4.24 

•  » 

Table  26.  Root  mean  square  error  (m)  for  48-hour  forecast  height  fields  valid 
at  00Z,  26  August  1983  for  the  control  run. 

Global  NH  SH 


Level  (mb) 

Obs. 

RMSE 

Obs. 

RMSE 

Obs. 

RMSE 

850 

665 

31.1 

588 

30.3 

77 

37.1 

500 

1412 

41.6 

1008 

34.8 

404 

55.1 

300 

1389 

55.9 

988 

49.5 

401 

69.3 

100 

1306 

75.5 

931 

70.5 

375 

86.6 

Table  27.  Root  mean  square  error  (m)  for  48-hour  forecast  height  fields  valid 
at  00Z,  26  August  1983  for  the  assimilation  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NH 

Obs.  RMSE 


SH 

Obs.  RMSE 


850 

658 

3.82 

581 

3.81 

77 

3.89 

500 

1416 

2.39 

1006 

2.26 

410 

2.69 

300 

1400 

2.49 

990 

2.41 

410 

2.68 

100 

1149 

5.36 

795 

5.45 

354 

5.15 

Table  29.  Root  mean  square  error  (°C)  for  48-hour  forecast  temperature  fields 
valid  at  00Z,  26  August  1983  for  the  assimilation  run. 


Level  (mb) 


Global 

Obs.  RMSE 


NH 

Obs.  RMSE 


SH 

Obs .  RMSE 


rror  (m)  for  72-hour  forecast  height  fields  valid 
1983  for  the  control  run. 


657 

34 

250 

51 

206 

69 

115 

82 

square  e 
7  August 


Table  32.  Root  mean  square  error  (°C)  for  72-hour  forecast  temperature  fields 
valid  at  00Z,  27  August  1983  for  the  control  run. 


Level  (mb 


Global 

Obs.  RMSE 


NH 

Obs.  RMSE 


SH 

Obs.  R 


Table  33.  Root  mean  square  error  (°C)  for  72-hour  forecast  temperature  fields 
valid  at  00Z,  27  August  1983  for  the  assimilation  run. 


Level  (mb 


Global 

Obs.  RMSE 


NH 

Obs .  RMSE 


646 

3.76 

573 

3.72 

73 

249 

2.84 

887 

2.67 

362 

223 

2.70 

859 

2.61 

364 

121 

3.86 

790 

3.60 

331 

forecasts  due  to  assimilation  is.  seen  at  all  levels  in  the  Northern  Hemi¬ 
sphere  but  only  at  100  mb  in  the  Southern  Hemisphere.  Again,  the  temperature 
forecasts  reflect  more  consistent  improvement  with  the  largest  seen  at  100  mb 
in  the  Northern  Hemisphere. 

In  summary,  the  results  from  the  second  experiment  were  a  little  dif¬ 
ferent  than  those  from  the  first.  The  improvements  in  100  mb  height  fore¬ 
casts  were  not  nearly  as  pronounced  for  the  second  experiment  although  the 
temperature  forecast  improvements  were.  Again,  consistent  improvement  was 
noted  in  the  temperature  forecasts  although  the  300  mb  level  again  showed 
little  improvement  or  degradation.  For  the  second  experiment  more  improve¬ 
ment  was  evidenced  in  the  Northern  Hemisphere  while  in  the  first  experiment 
the  largest  improvements  were  seen  in  the  Southern  Hemisphere.  We  can  only 
speculate  about  the  possible  causes  for  these  differences.  One  reason  mig>"t 
be  the  difference  in  seasons  as  the  first  experiment  was  run  in  mid-May  while 
the  second  was  run  in  late  August.  The  rmse  calculations  were  made  using 
all  available  observations  (conventional  and  satellite)  but  that  still  does 
not  provide  truly  global  coverage.  There  very  well  could  be  fluctuations  in 
these  rmse  calculations  due  to  the  extent  of  the  observational  coverage.  In 
any  case,  however,  it  is  encouraging  to  note  that  in  both  experiments  the 
initial  fields  from  the  assimilation  runs  showed  rms  improvement  over  those 
from  the  control  runs  and  that  temperature  forecasts  also  showed  consistent 
and  often  marked  improvement.  The  fact  that  the  height  fields  did  not  always 
display  the  same  improvement  that  the  temperature  fields  did  indicates  that 
some  attention  needs  to  be  paid  to  the  pressure  forecasts  made  during  the 
assimilation  runs. 

The  forecast  fields  for  24,  48,  and  72  hours  from  00Z,  24  August  1983 
are  displayed  in  Figs.  3-62  for  both  the  control  run  and  assimilation  run. 
Comparing  Figs.  3  and  13,  23  and  33,  and  43  and  53  we  see  that  only  subtle 


differences  between  the  two  runs  are  apparent  for  forecasts  of  sea-level 
pressure  in  the  Northern  Hemisphere.  In  Figs.  4  and  14  we  see  that  after 
24  hours  the  assimilation  forecast  for  sea-level  pressure  depicts  lows  at 
about  60S,  90E  and  at  about  50S,  SLOW  to  be  deeper  than  the  control  forecast. 
Also  the  high  latitude  low  at  around  150W  is  portrayed  quite  differently. 

After  48  hours  we  see  in  Figs.  24  and  34  that  the  low  at  60S,  90E  is  only 
slightly  deeper  now  but  that  the  low  off  the  tip  of  South  America  is  deeper 
and  defined  differently.  Again  the  high  latitude  low  near  150W  is  much 
different  and  a  high  has  developed  in  the  control  forecast  at  about  40S,  110W. 
Finally,  an  inverted  trough  at  about  30S,  180W  is  more  pronounced  for  the 
assimilation  run.  In  Fig.  44  and  54  we  see  that  the  lows  at  60S,  90E  are 
defined  quite  differently  as  is  the  high  latitude  low  near  150W  amd  the  low 
off  the  tip  of  South  America.  The  high  near  40S,  105W  is  much  stronger  for 
the  control  run  while  the  assimilation  run  has  developed  a  deeper  low  near 
30S,  170W.  Thus,  we  have  seen  in  these  figures  that  especially  in  the 
Southern  Hemisphere  significant  qualitative  differences  in  the  depiction  of 
meteorological  features  have  developed  between  the  control  and  assimilation 
forecasts  for  sea-level  pressure.  In  this  study  we  do  not  have  the  information 
necessary  to  assess  the  quality  of  these  forecasts  but  it  is  important  to  note 
the  differences  that  occur.  Only  by  daily  operational  testing  could  one 
expect  to  determine  which  run  (control  or  assimilation)  actually  produces 
better  forecasts  of  these  meteorological  features. 

Looking  at  the  forecast  of  850  mb  heights  in  the  Northern  Hemisphere,  we 
see  in  Figs.  5  and  15  that  after  24  hours  the  assimilation  run  portrays  the 
structure  of  the  low  near  Kamchatka  differently  along  with  that  of  the  trough 
off  the  west  coast  of  the  U.S.  A  low  in  the  area  of  the  Alps  is  also  deeper 


for  the  assimilation  run.  After  48  hours  we  can  still  see  differences  in 
Figs.  25  and  35  between  these  features  for  the  two  runs  although  they  appear 
to  be  more  subtle.  In  Figs.  45  and  55  we  see  that  after  72  hours  the  low 
off  Kamchatka  is  depicted  by  the  assimilation  run  to  be  deeper  with  a  dis¬ 
tinct  trough  approaching  the  dateline,  a  low  center  has  developed  in  the 
Alaskan  Panhandle  area,  and  only  minor  differences  are  evident  in  the  hand¬ 
ling  of  the  Alpine  low.  In  the  Southern  Hemisphere  Figs.  6  and  16  show 
that  after  24  hours  differences  in  structure  are  seen  between  the  two  runs 
in  the  handling  of  the  low  at  about  60S,  90 E,  the  high  southwest  of  New 
Zealand,  the  troughing  around  150W,  and  the  low  off  the  western  tip  of 
South  America.  After  48  hours  we  see  in  Figs.  26  and  36  that  the  most 
notable  differences  are  in  the  forecasts  for  the  trough  around  150W  and  the 
low  off  the  tip  of  South  America.  For  the  72-hour  forecasts  we  see  in  Figs. 
46  and  56  that  the  assimilation  run  develops  a  relatively  deep  low  center 
at  about  30S,  170W  while  the  control  run  displays  only  weak  troughing. 

The  trough  at  about  140W  and  the  low  off  the  tip  of  South  America  are  also 
portrayed  quite  differently. 

As  we  look  higher  in  the  atmosphere  the  differences  between  the  assimi¬ 
lation  and  control  forecast  fields  become  more  and  more  subtle.  In  Figs. 

7  and  17  the  only  really  noticeable  difference  between  the  24-hour  Northern 
Hemisphere  500  mb  forecasts  is  in  the  depiction  of  the  trough  off  the  west 
coast  of  the  U.S.  and  the  low  near  Kamchatka.  We  can  make  the  same  observa¬ 
tions  for  the  48-  and  72-hour  forecasts  displayed  in  Figs.  27  and  37  and 
Figs.  47  and  57,  respectively.  The  same  can  be  said  for  the  Southern  Hemi¬ 
sphere  where  we  see  in  Figs.  8,  18,  28,  38,  48,  and  58  that  the  only  notice¬ 
able  differences  are  in  the  treatment  of  the  trough  east  of  New  Zealand,  the 


trough  around  150W  and  the  trough  off  the  tip  of  South  America.  In  the  re¬ 
maining  figures  we  see  that  for  the  forecasts  at  30Q  mb  and  100  mb  only  subtl 
differences  are  apparent  with  a  few  exceptions.  Comparing  Figs.  10  and  20 
we  see  that  the  trough  off  the  South  American  tip  at  300  mb  is  much  sharper 
for  the  24-hour  assimilation  forecast  and  the  low  at  about  60S,  90E  is 
closed  off.  In  Figs.  30  and  40  the  trough  off  the  tip  of  South  America  is 
still  sharper  for  the  48-hour  assimilation  forecast  and  the  same  holds  true 
after  72-hours  as  evidenced  by  Figs.  50  and  60. 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  the  previous  section  we  described  the  results  of  two  continuous  data 
assimilation  experiments  conducted  using  NOGAPS  and  the  satellite  data  opera¬ 
tionally  available  to  FNOC.  The  techniques  used  in  these  experiments,  des¬ 
cribed  in  Section  2,  were  tailored  for  application  with  NOGAPS.  In  this 
section  we  shall  summarize  our  findings  and  make  recommendations  for  the 
future  application  of  our  research. 

The  first  important  result  we  obtained  is  that  the  NOGAPS  model  is 
ideal  for  continuous  data  assimilation.  Over  the  course  of  our  data  assimi¬ 
lation  research  we  have  found  that  one  way  to  obtain  forecast  improvement  is 
to  control  the  model  noise  induced  by  data  insertion.  The  Noise  Freezing 
Methods  described  in  Sasaki  and  Goerss  (1982)  were  developed  to  accomplish 
this.  We  saw  in  Section  2  that  the  Matsuno  time-step  which  initiates  each 
integration  cycle  of  the  NOGAPS  model  has  a  damping  effect  upon  this  noise. 
Thus,  we  found  that  we  were  able  to  continuously  assimilate  temperature 
data  into  the  model  without  having  to  make  any  changes  to  the  model  itself. 

The  variational  adjustment  scheme  we  call  vertical  shape  matching  was 
designed  to  preserve  the  information  content  of  each  satellite  sounding  and 
yet  modify  it  to  make  it  easier  for  the  forecast  model  to  digest.  Its 
success  was  demonstrated  in  Section  2  where  we  saw  that  the  damping  prop¬ 
erties  of  the  NOGAPS  model  kept  the  noise  under  control  while  the  temperature 
data  was  assimilated  into  the  model  without  any  indications  of  rejection. 

We  described  in  Section  3  one  of  the  most  striking  results  of  continu¬ 
ous  assimilation.  The  global  mean  temperatures  for  the  upper  two  layers  of 
the  NOGAPS  model  showed  marked  changes  due  to  assimilation  during  the  start-up 
period  for  each  experiment.  After  the  two-day  start-up  period,  the  upper 


level  field  for  the  assimilation  run  was  about  10  C  colder  than  that  for  the 
control  run.  Recall  that  the  satellite  data  was  also  used  in  the  control  run 
but  only  in  the  objective  analysis  scheme  applied  each  update  cycle  (every 
six  hours).  At  the  next  lower  level  the  mean  temperature  for  the  assimila¬ 
tion  run  was  about  3°C  warmer  than  that  for  the  control  run.  This  modifica¬ 
tion  to  the  mean  NOGAPS  model  sounding  caused  by  continuous  data  assimilation 
resulted  in  a  more  accurate  depiction  of  the  conditions  in  the  actual  atmos¬ 
phere. 

Next,  forecasts  made  from  the  assimilation  and  control  runs  for  each 
experiment  were  quantitatively  evaluated  by  computing  root  mean  square  errors 
for  the  differences  between  actual  observations  of  geopotential  height  and 
temperature  and  the  values  of  forecast  fields  interpolated  to  the  location 
of  the  observations.  In  general,  improvement  was  noted  due  to  assimilation. 
The  most  significant  improvements  were  seen  at  100  mb  for  both  geopotential 
height  and  temperature.  The  improvements  for  temperature  were  much  more 
widespread  and  marked  than  those  for  height.  Differences  were  noted  for  the 
results  between  the  two  experiments  which  could  have  been  seasonal  or  due  to 
the  observational  coverage.  In  summary,  however,  the  results  of  the  rmse 
computations  implied  that  continuous  data  assimilation  had  a  positive  impact 
overall  on  the  NOGAPS  forecasts. 

Finally,  a  qualitative  assessment  of  the  forecasts  was  made  for  the 
second  experiment.  Conclusive  statements  could  only  be  made  about  these 
results  if  the  assimilation  and  control  runs  were  compared  on  a  daily  basis 
in  an  operational  evaluation  over  a  significant  period  of  time.  However,  in 
this  study  we  found  that  noticeable  differences  were  seen  in  the  handling 
of  meteorological  features  by  the  two  runs,  especially  in  the  Southern 
Hemisphere. 


From  our  research  we  conclude  that  the  forecasts  produced  by  NOGAPS  can 
be  improved  by  the  use  of  continuous,  data  assimilation.  The  operational 
expense  of  this  technique  would  be  that  required  to  run  an  additional  six- 
hour  assimilation  forecast  each  update  cycle.  It  was  very  encouraging  to 
discover  that  the  NOGAPS  model  lends  itself  so  well  to  continuous  assimilation 
without  modification.  There  are  very  few  parameters  that  need  to  be  consid¬ 
ered  in  order  to  fine  tune  the  continuous  data  assimilation  techniques 
described  in  Section  2  for  application  in  a  truly  operational  setting. 

Minor  adjustments  to  the  weighting  parameters  in  the  indirect  insertion 
scheme  (which  was  patterned  after  that  described  by  Barker  (1982)  and  is 
similar  to  the  objective  analysis  scheme  used  in  NOGAPS)  could  be  made  and 
forecasts  evaluated.  Special  attention  needs  to  be  paid  to  the  effects  of 
assimilation  upon  the  forecasts  of  surface  pressure  during  this  "tuning" 
process.  Although  the  impact  of  continuous  satellite  data  assimilation 
upon  the  NOGAPS  forecasts  can  only  be  determined  after  extensive  operational 
testing,  the  results  of  our  research  reported  here  indicate  that  significantly 
positive  impact  can  be  expected. 
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FORECAST  TIME  (hr) 


(Top)  Root  mean  square  difference  (°C)  between  satellite 
temperature  observations  at  592  mb  and  NOGAPS  model 
temperatures. 

(Bottom)  Global  root  mean  square  pressure  tendencies 
(mb/sec)  for  assimilation  (top)  and  non-assimilation 
NOGAPS  forecasts. 


Fig.  10.  24-hour  Southern  Hemisphere  assii 
heights  (m)  valid  at  00Z,  25  Aug 


24-hour  Northern  Hemisphere  control  run  forecast  for  sea-level 
pressure  (mb)  valid  at  00Z,  25  August  1983. 


24-hour  Northern  Hemisphere  control  run  forecast  for  500  mb 
heights  (m)  valid  at  00Z,  25  August  1983. 


24-hour  Northern  Hemisphere  control  run  forecast  for  100  mb 
heights  (m)  valid  at  00Z,  25  August  1983. 
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Fig.  26.  48-hour  Southern  Hemisphere  assimilation  run  forecast  for  850  mb 
heights  (m)  valid  at  00Z,  26  August  1983. 
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28.  48-hour  Southern  Hemisphere  assimilation  run  forecast 
heights  (m)  valid  at  00Z,  26  August  1983. 


hour  Southern  Hemisphere  assimilation  run  forecast  for  300  mb 
ghts  (m)  valid  at  00Z,  26  August  1983. 


Fig.  32.  48-hour  Southern  Hemisphere  assimi 
heights  (m)  valid  at  00Z,  26  Augus 


Fig.  36.  48-hour  Southern  Hemisphere  control  run  forecast  for  850  mb 
heights  (m)  valid  at  00Z,  26  August  1983. 
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27  August  1983 


72-hour  Southern  Hemisphere  assimilation  run  forecast  for  500  mb 
heights  (m)  valid  at  00Z,  27  August  1983. 
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72-hour  Southern  Hemisphere  assimilation  run  forecast  for  300  mb 
heights  (m)  valid  at  00Z,  27  August  1983. 
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72-hour  Northern  Hemisphere  assimilation  run  forecast  for  100  mb 
heights  (m)  valid  at  00Z,  27  August  1983. 


83 

.v.v.v 


Fig.  52.  72-hour  Southern  Hemisphere  assimi 
heights  (m)  valid  at  00Z,  27  Augus 
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72-hour  Northern  Hemisphere  control  run  forecast  for  sea-level 
pressure  (mb)  valid  at  00Z,  27  August  1983. 


Fig.  58.  72-hour  Southern  Hemisphere  cont 
heights  (m)  valid  at  00Z,  27  Augi 
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